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Abstract 

In this article we propose a simultaneous collider search strategy for a pair of scalar 
bosons in the NMSSM through the decays of a very light pseudoscalar. The massive 
scalar has a mass around 126 GeV while the lighter one can have a mass in the 
vicinity of 98 GeV (thus explaining an apparent LEP excess) or be much lighter. 
The successive decay of this scalar pair into two light pseudoscalars, followed by 
leptonic pseudoscalar decays, produces clean multi-lepton final states with small or 
no missing energy. Furthermore, this analysis offers an alternate leptonic probe 
for the 126 GeV scalar that can be comparable with the ZZ* search channel. We 
emphasize that a dedicated experimental search for multi-lepton final states can 
be a useful probe for this scenario and, in general, for the NMSSM Higgs sector. 
We illustrate our analysis with two representative benchmark points and show how 
the LHC configuration with 8 TeV center-of-mass energy and 25 fb _1 of integrated 
luminosity can start testing this scenario, providing a good determination of the 
light pseudoscalar mass and a relatively good estimation of the lightest scalar mass. 
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1 Introduction 



The experimental evidence for a scalar boson with a mass around 126 GeV [H 
[2] has been received with enthusiasm by the particle physics community, since it 
could correspond to the missing piece of the standard model (SM), the Higgs boson. 
Although most of the properties of this new particle (in particular the branching 
ratios of various decay modes) are compatible with a SM-like Higgs, an apparent di- 
photon excess [3|H] . has been considered as a possible hint for physics beyond the SM, 
such as supersymmetry (SUSY). However, this excess remains yet to be confirmed [5]. 
It has been pointed out that in order to reproduce the mass of the Higgs boson in 
the Minimal Supersymmetric Standard Model (MSSM) one generally needs a very 
heavy spectrum, raising concerns about the naturalness of this construction. 

This problem can be alleviated in extended Supersymmetric constructions. In 
particular, the Next-to-Minimal Supersymmetric Standard Model (NMSSM) can 
accommodate a 126 GeV scalar [BTfTT]. while maintaining the idea of naturalness 
[IlT[l~4"] . The advantage of the NMSSM resides in the presence of an extra singlet 
field, which provides an additional contribution to the tree- level Higgs mass [ToTflT] . 
Thus, contrary to what happens in the MSSM, loop corrections can be smaller and 
the correct Higgs mass is achieved for a lighter supersymmetric spectrum. This 
extra singlet superfield S is introduced to promote the bilinear /i-term nHdH u in 
the MSSM to a trilinear coupling \SHdH u . After electroweak symmetry breaking 
(EWSB) takes place the vacuum expectation value (VEV) of the singlet field triggers 
an effective [i parameter \v s which is naturally of order of the electroweak scale, thus 
curing the //-problem of the MSSM [18]. On top of this, the singlet-doublet mixing of 
the resulting scalar states is useful to produce the observed di-photon excess [6lfT9j. 

The new singlet superfield induces one extra CP-even and CP-odd neutral scalars 
as well as one additional neutral fermion, thereby leading to a very rich phenomenol- 
ogy. In particular, a light pseudoscalar Higgs is viable if its singlet component is 
large [20ti27] . Similarly, the lightest supersymmetric particle (LSP) can also be 
singlino-like [22] or, in general, a mixed state with interesting implications for dark 
matter searches (see, e.g., Refs. [29||30]). 

The presence of light scalar and pseudoscalar Higges can lead to characteristic 
Higgs-to-Higgs cascade decays [20j[2T]. Although these cascades are generally sup- 
pressed, a substantial enhancement can occur when particles are light and therefore 
induce distinctive signatured^- For example, if the lightest pseudoscalar, denoted 
as a 1 ] 1 , has a mass m a o > 2mb, the process h\ — > 2a® — > 2bb with 46 final states is 

1 To some extent, a scenario like this is a good illustration of the "no-lose" theorem for the 

nmssm H0H21I3IH3S]. 
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dominant [2"2" |l3"6T - r4"0"] . This seems a good channel for collider studjo, however, in 
hadron colliders like the LHC this signature is hidden by a huge QCD background. 
A viable alternative is looking for At [201 [221 [251 EU [351 [391 SU S2] or 2/i2r [42] final 
states, for which the problem of large QCD backgrounds is somehow ameliorated. 
However, this approach suffers from the poor tau-identification efficiency, dependent 
on the transverse momentum [43]. Besides, since tau decays are always accompanied 
by neutrinos, there is no sharp peak of a? invariant mass. For these reasons, as dis- 
cussed in the literature, 4/i [39HH1I15] , Aj [2T1[33] or 27 + 2-jets [46] search channels 
with m a o < 2m T turn out to be the ideal way to probe NMSSM Higgses with very 
light pseudoscalars. 

Recently, the possibility of having two light scalar Higgses in NMSSM has re- 
ceived a lot of attention. One of these Higgses, h®, would correspond to the scalar 
observed by ATLAS and CMS with a mass around 126 GeV and the other one, h®, 
with a mass around 98 GeV pU], would be consistent with the small excess in the 
LEP search for e + e~ — > Zh, h — >■ bb [47] |§ Alternative scenarios with an even lighter 
hi have also been considered earlier [22l[50]. Since the recent LHC observation of 
the 126 GeV scalar, which is well compatible with a SM-like Higgs boson, many 
of these studies have become extremely constrained, although a small window for 
new physics is still open if the di-photon excess is confirmed. On the other hand, 
scenarios with very light pseudoscalars are also affected by stringent experimental 
bounds. This is the case for m a o < 2m T (which leads to leptonic final states 4/i, 
2/i2t, and At) [261 [Ml 021 [51] > an d also for m a o < 2m^ which is tightly constrained 
for CP-even Higgs masses in the range of 86 - 150 GeV by the recent CMS analy- 
sis [52], especially in the case of singlet-like states. The latter is a consequence of 
the fact that BR(hl/hl —> a^a?) ~ 1 and thus a(pp -> h\jh\ — > 2a\ — >■ 4/i) is sizable 
and would have already been observed (a way to avoid this is if the production cross 
section for a singlet-like h\jh\ is very suppressed). 

In this work we investigate potential detection channels for scenarios involving 
very light pseudoscalar particles and two light scalar Higgses in the NMSSM. For 
concreteness, we consider a SM-like scalar Higgs in the range 124 GeV < m h a < 128 
GeV (consistent with LHC findings) and a lighter h\ in the mass range 96 - 100 GeV 
(to account for the LEP excess) or lighter than 86 GeV (to avoid the current CMS 
limit). In order to avoid the above mentioned constraints that affect the lightest 
pseudoscalar, a\ is assumed to have a mass in the range 2m T < m a o < 2m&, since 
the final state is then a sub-leading leptonic mode with a branching ratio 

of the order of 10~ 2 and a(pp — > h\ — > 2a\ — > 4/i) is small enough to elude the 
CMS limit of the 4/i search. We point out that the SM-like Higgs decay modes 

2 Another possible search channel is 262r. 

3 A two light Higgs but without a low-mass pseudoscalar has been discussed in the MSSM [48] . 
as well as in the NMSSM gS]. 
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h® — > 2a5 — > 4/i, 2/i2r, and 4r have branching fractions comparable to those for 
the same final states from /t° — * ZZ* and therefore provide alternative four-lepton 
signatures (accompanied with zero or a small missing energy) at the LHC. We also 
note that the long Higgs-to-Higgs cascades of h® — > h\h\ — > Aa\ can be open if 
m h o < 86 GeV, which constitutes a clean and distinctive multi-lepton signal at the 
LHC. For these channels, we study the distributions of the kinematic variable M^ 2 
and the invariant mass of the four visible leptons, and determine to what extent 
the masses of the pseudoscalar and scalar Higgs particles can be reconstructed. In 
our analysis we include constraints from i?-physics [531456] . the muon anomalous 
magnetic moment [5TH60] . and the LHC result for h% — > 77. In addition, we impose 
the WMAP bound on the lightest neutralino relic abundance [61] and we ensure that 
its spin-dependent scattering cross section is consistent with XENON100 data [62] . 

The paper is organized as follows. We start with a brief description of the model 
and determine the relevant signatures in Sec. [2] The choice of model parameters is 
justified in Sec.[3]where we also describe the details of the analysis. Sec.H]is dedicated 
to the analysis of the numerical results. We also introduce some discriminating 
variables to suppress SM backgrounds and we investigate the statistical significance 
of the proposed signals. Finally we summarize our conclusions in Sec. [51 



2 Higgs signals with a low-mass pseudoscalar in 
the NMSSM 

The superpotential of the Z 3 -symmetric NMSSM is given by 

W — W — e ab \SH a d H b u + ^kSSS, (1) 

where W is the MSSM superpotential without the /i-term, S is a superfield singlet 
under the SM gauge group, and a, b = 1, 2 are SU(2) L indices. Similarly, with £g oft 
representing the MSSM soft-terms without the -B M -term, the Lagrangian containing 
the soft SUSY-breaking terms in a supergravity framework is given by 

- £ soft = -4 0ft + m 2 s S*S - e ab (\A x )SH a d H b u + 1 -(kA k )SSS + H.c. . (2) 

The last trilinear term in Eq. ([1]) is essential to avoid an unacceptable axion asso- 
ciated with the breaking of a global U(l) symmetry [16]. After EWSB, this term 
also generates an effective Majorana mass 2kv s for the singlino-like neutralino. The 
second term of Eq. flTJ not only generates an effective \i parameter, /i eff = Xv s , but 
also provides an extra tree-level contribution to the lightest doublet-like Higgs boson 
mass. The complete expression is now [T5l - [T7] . 

m 2 h0 < M 2 Z cos 2 2(3 + 3.62 Aff A 2 sin 2 2(3, (3) 
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where tan/3 = v u /v d is the ratio of up and down- type Higgs VEVs. Thus even 
maintaining perturbativity (A < 0.7), the extra contribution to the tree-level Higgs 
mass can be sizable for small values of tan /?. The last term of Eq. ([2]) is important 
to understand the phenomenology of singlet-like scalars. In particular, when the 
lightest CP-even scalar h\ as well as the lightest CP-odd scalar a\ are predominantly 
singlet-like, the decay width for h® — > a\a\ is proportional to A k k. 

As already mentioned in the previous section, the Higgs sector of the NMSSM is 
richer than that of the MSSM, featuring three CP-even scalar states, (h®, h®, h®) and 
two pseudoscalars, (a°, a^). This is of particular interest for studying Higgs cascade 
decays, which can now be long and often result in multi-particle final states with no 
missing energy. The decay chains become more interesting in the presence of light 
scalar and pseudoscalar Higgses, which can be in agreement with present constraints 
if their singlet composition is dominant. In particular, a viable hierarchy for the 
scalar mass eigenstates in the NMSSM consists of having a light singlet-like h®, a 
SM-like Higgs h® with a mass around 126 GeV and a much heavier h® which is also 
doublet-like. This structure with two light scalar Higgses might produce distinctive 
features in collider searches, as compared to the signals that come from the MSSM 
Higgs sector, and therefore we concentrate on this possibility. 

On top of this, very light singlet-like pseudoscalars can give rise to potentially 
characteristic signatures, but one must be very careful with recent experimental 
constraints. In particular, for a scalar Higgs mass in the mass range of 86 - 150 GeV, 
the upper limit of a(pp —> h^/h^ — >■ 2a®) xBR 2 (a® —> has been found to be 3 

- 5 fb, depending on the pseudoscalar mass for 1 GeV < m a o < 2m T in the 7 TeV 
CMS search [52]. Thus, if we want to study the 98 + 126 GeV scenario for the scalar 
Higgs bosons, in order to avoid this constraint, we will consider a light pseudoscalar 
with a mass in the range 2m r < m a o < 2m b . We call this Scenario I. 

Another way to avoid the CMS bound is to consider a lighter scalar Higgs with 
m h o < 86 GeV. Notice that although this would relax the LHC limit on the pseu- 
doscalar mass, in fact m a o < 2m T is still severely constrained from the LEP limit 
on 4r or multi-lepton/jet final states |36j|63] in the case when m h o < m h o/2. Thus, 
2m T < m a o < 2m b is still the preferred range for pseudoscalar masses. We call this 
Scenario II. 

For clearness we summarise here the properties of the two kind of scenarios 
considered in our work 

(I) 124 GeV< m h o <128 GeV, 
96 GeV< m h o <100 GeV, 
2m T < m a o < 2m b , 

(II) 124 GeV< m h o <128 GeV, 
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m h o < m h o/2, 
2m T < m a o < 2m b . 

We have also checked that these benchmark scenarios are in agreement with ATLAS 
and CMS direct pseudoscalar searches, pp — > a° — > fi + fi~ [M], mainly due to the 
singlet nature of a\. 

We are now ready to specify the potential signatures that we will study for these 
scenarios. We consider the following decay modes, 

/i°^a°a°-> 4£ + ^ T , (4) 

where I = e, /i, r. We here assume that h\ and h\ have been produced through 
gluon-fusion and consider inclusive decay modes of the tau lepton. The missing 
energy is associated with neutrinos coming from tau-decays and is generally 
small. It is important to emphasize that direct-pseudoscalar decays to a muon-pair 
will eventually yield a clean four-muon final state with vanishing missing energy and 
the only sources of electrons are leptonic r-decays, since BR(a5 —> e + e~) ~ 0. The 
second decay mode can act as an alternative source of four-lepton final state apart 
from h® — > ZZ* — > 4£. In the chosen corner of parameter space, we found that the 
predicted BR(/i2 — > a? a?) is comparable with that of BR(/i2 — > ZZ*), which still 
matches well the SM prediction. We thus encourage our experimental colleagues 
to search for possible excesses in 4/x, 4r or 2/i2r channels by relaxing the Z-boson 
resonance condition. In the remainder of the paper we carry out a dedicated analysis 
including the discussion of SM backgrounds and useful collider variables to search 
for the four-lepton collider signals of these decay modes. 

Concerning possible SUSY backgrounds to these decay modes, there can be many 
sources of four-lepton final states with a small or no missing energy, provided that 
one considers -R-parity violation in the MSSM [65], NMSSM [66], or /ii^SSM [67] . 
It is well known that -R-parity conserving SUSY models are characterized by large 
missing energy due to the production of a heavy neutralino LSP, which is stable and 
escapes the detector, whereas this is not the case in i?-parity violating scenarios. On 
the other hand, all these models exhibit moderate or large displaced vertices with 
non-prompt jets or leptons in the final states (see e.g., Refs. [SHIEH]), more complex 
Higgs cascade decays [69], or occasional correlations with neutrino physics [70] that 
can be used to distinguish among them. We also note that four leptons can be 
produced in ft^ decays, but with a very reduced production cross section. 

Notice that one distinguishing feature of Scenario II is a three-step cascade decay 
of the Higgs boson h\ — > 2h\ — > Aa\ — > multi-lepton/jet. We have checked for several 
points in the parameter space that, even when m h o > m h o/2, the branching ratio of 
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h® — > — > Aa\ can be as large as ~ 10 -3 and thus multi-lepton final states 
can be experimentally accessible. Due to the large number of leptons and/or jets 
in the final state, the main backgrounds are expected to be SUSY processes (such 
as cascade decays mediated by neutralinos and/or charginos), rather than SM ones. 
The analysis of this decay signal will be carried out in a separate work [71] since the 
search strategy is very different than for the other two. 



3 Experimental constraints and choice of param- 
eters 

The recent discovery of a 126 GeV scalar along with the reported di-photon excess 
set stringent limitations on the NMSSM parameter space. In addition, considering 
bounds from low-energy observables in flavour physics, the supersymmetric contri- 
bution to the muon anomalous magnetic moment a^ USY , and the relic abundance 
of the dark matter, which in our case would correspond to the lightest neutralino, 
narrow down significantly the viable regions. 

In our analysis we impose compatibility with all these constraints. The numer- 
ical results are obtained with nmssmtools 3.2.1 [72], which calculates the mass 
spectrum and provides predictions for low-energy observables, as well as computing 
the decay widths of Higgses [73] and sparticles [73] . The masses of Higgs bosons 
include full two-loop contributions. We also consider the lightest neutralino as a 
dark matter candidate and include the WMAP upper constraint on its relic abun- 
dance [61] as well as the upper bound on the spin-independent neutralino-nucleon 
cross section from XENON100 [H2]- The neutralino relic density and its scattering 
cross section are computed through an interface with micromegas [75]. We have 
also modified the code nmssmtools to include three-body Higgs decay in the case 
of the NMSSM. 

We have calculated the theoretical predictions for BR(£>° — > and com- 

pared them with the experimental value. In the case of the MSSM, the SUSY 
contributions to this observable [76] are likely to exceed the recent LHCb measure- 
ment BR(B S ° — > /i + yU~)= 3.2 x 10~ 9 [55], especially for light pseudoscalars in the 
large tan/3 regime. In the case of the NMSSM [77J,[78], the corresponding Wilson 
coefficient receives contributions from both pseudoscalar Higgs bosons through their 
doublet components. Given that the light pseudoscalar is a purely singlet-like field, 
it leads to a negligible effect on this observable, and since a° is heavy enough and 
tan /3 is small, the experimental constraint is easily fulfilled. On the other hand, in 
general, BR(6 — > sj) provides a strong constraint on SUSY models. On top of the 
usual MSSM terms, NMSSM-specific contributions arise from the extended Higgs 
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Parameter 



BP1 



BP2 



tan /3 
A, K 



0.285, 0.1165 
670, 14.0 
300, 300 



5 



0.286, 0.0844 
820, 14.35 



5 



A x , A K (GeV) 
M z ., M g c (GeV) 
Mq/m^Mq (GeV) 



1000, 1000, 1000 

123.5 
560, 1200, 1980 
-1600, 1000, 1800 



300, 300 
1000, 1000, 1000 



IM (GeV) 



123.5 
240, 500, 1380 
-1600, 1000, 1300 



Mi, M 2 , M 3 (GeV) 
A,, A b , A t (GeV) 



Table 1: Model parameters that define our choice of benchmark points. The top-quark 
pole mass is set to 173.5 GeV and mt, us (mi,) = 4.18 GeV. 

and the neutralino sectors, which come into effect at the two-loop level [77]. This 
observable has been shown to lead to stringent constraints on the NMSSM parame- 
ter space [30] and this is indeed the case in our current analysis. Here, we consider 
the experimental result [56] and include the theoretical error of the calculation in 
the SM [79] in quadrature. 

The SUSY contribution to the muon anomalous magnetic moment, a^ USY was 
also computed. The observed discrepancy between the experimental value [57] and 
the SM predictions using e + e~ data favours positive contributions from new physics 
in the range 10.1 x 10~ 10 < c^ t USY < 42.1 x 10" 10 at the 2a confidence level [581 - 
[60], combining experimental and theoretical errors in quadrature. If tau data is 
employed, this discrepancy is smaller, 2.9 x 10~ 10 < a^ USY < 36.1 x 10~ 10 [59] . 

We have carried out a simple scan to sample the NMSSM parameter scan, search- 
ing for regions in the parameter space where Scenarios I and II could be realised. 
Out of the viable regions we have selected two benchmark points, BP1 and BP2 
with parameters given in Table [TJ which constitute representative examples of Sce- 
nario I and Scenario II, respectively. The first two generation squarks and gluino 
masses are chosen to be heavy enough to be in agreement with current LHC SUSY 
searches [80]. We have relaxed gaugino mass unification in order to have more 
freedom in the neutralino composition, which is important in order to fix its relic 
abundance and scattering cross-section. Also, trilinear parameters are non-universal 
and chosen in such a way that a^ USY is maximized and BR(6 — > 37) is in agreement 
with the experimental value. The resulting mass spectra are shown in Table [3 

The compositions of Higgses, neutralinos, and charginos for our benchmark points 
are shown in Table O Since the \x parameter is chosen to be small and A ~ 0.3, the 
singlet VEV is around 430 GeV. Thus with the small k value, the singlino mass, 
2kv s for the benchmark points BP1 and BP2 is approximately 72 and 100 GeV, 
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Parameter 


BP1 


BP2 


m h o, m h o, m h o 


97.7, 125.9, 677.0 


62.1, 125.9, 739.8 


"la". m a% 


3.6, 675.4 


7.6, 738.6 


m h ±, rn~±, m~± 


679.8, 124.0, 1192 


739, 118.1, 522.3 


m^o. m^o, m^o 

Xi' A 2 ' A3 


87.0, 136.1, 143.3 


63.6, 124.8, 139.0 


m ~t^ ™V m ~ b2 


859.5, 1154, 1008, 1009 


951.7, 1139, 1040, 1041 




296.5, 309.5 


296.5, 309.4 


rrig 


1926 


1393 



Table 2: Relevant mass spectrum for the chosen benchmark points. For h® the latest 
ATLAS limit is 125.2 ± 0.3(stat) ± 0.6(syst) GeV [81]. The CMS limit corresponds to 
125.3 ± 0.4(stat) ± 0.5(syst) [2]. Masses are given in GeV. 



Mass eigenstate 


BP1 (%) 


BP2 (%) 


hi (S, H u , H d ) 


79.9, 17.6, 2.5 


81.2, 16.5, 2.3 


h% (S, H u , H d ) 


19.5, 8,2, 2.3 


18.3, 79.4, 2.3 


h° 3 


95.2% H d -like 


95.4% H d -\ike 


a? 


99.6% singlet-like 


99.5% singlet-like 


a° 


95.8% ifd-like 


95.6% # d -like 


X\ (5°, W$, H d , H u , S) 


0.3, 0.2, 13.3, 30.0, 56.2 


1.7, 0.9, 7.8, 22.7, 66.9 


X2 (Higgsino, singlino) 


96.8, 2.9 


64.0, 29.5 


X3 (Higgsino, singlino) 


58.8, 40.9 


95.6, 3.6 


xt 


> 95% Higgsino-like 


> 95% Higgsino-like 



Table 3: Compositions of Higges, neutralinos, and charginos for the chosen benchmark 
points. 



respectively. Notice that this value is close to that of the fi parameter and that 
both of them are smaller than the bino and wino mass terms. This results in an 
interesting hierarchical spectrum, especially for the three lightest neutralino states. 
The lightest neutralino is a singlino-Higgsino state (where the singlino component 
slightly dominates). There is an orthogonal eigenstate which is also singlino-Higgsino 
(but this one with a larger Higgsino component), which in BP1 corresponds to the 
third neutralino and in BP2 is the second, and close in mass to this one we can find 
an almost pure Higgsino (x° f° r BP1 and \% f° r BP2). The heavier states are bino 
and wino- like. 

The resulting predictions for low-energy observables are shown in TableHJ Al- 
though the values of a^ USY are in the 2a range of the result from tau data, there is 
some tension with the one from e + e~ data. We have tried reducing this discrepancy 
by decreasing the masses in the slepton sector. Regarding the neutralino relic abun- 
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Quantity 


BP1 


BP2 


Range 


BR(6 s 7 ) x 10 4 


3.86 


3.69 


2.86 - 4.24 (2a) [56] 


BR(B S x 10 9 


3.97 


3.63 


2.0-4.7 [55] 


BR(B+ -»■ r+z/ T ) x 10 4 


1.31 


1.31 


0.85 - 2.89 (2a) [82] 


i? 77 


1.0295 


1.0385 


1.80 ±0.50 [ATLAS] [HE] 

1.4818;** [cms] m 


a SUSY x 1Q 10 


3.134 


6.90 


2.9 - 36.1 (2a) [59] 
10.1 - 42.1 (2a) [60] 


Relic density (Q^oh 2 ) 


0.082 


0.102 


0.094 - 0.136 [WMAP] [61] 


a SI x 10 9 (pb)' 


1.011 


0.723 


< 1.44@86.99 (2a) (BP1) [52] 

< 1.17@63.6 (2a) (BP2) [62] 



Table 4: Low-energy observables for the chosen benchmark point. When scanning the 
parameter space, we set -R 77 = J(^g^hsM-^ni) > ^ CSE]> within the 2a range of the 
ATLAS and the CMS result [IHlj. SI represents spin- independent. For BR(6 — > sj), the 
theoretical error is added to the experimental one in quadrature. 

dance, it is well known that a Higgsino-like neutralino neutralino generally entails a 
large annihilation cross section and consequently a small relic abundance, but in our 
case, the presence of a sizable (slightly dominant as we said above) singlino compo- 
nent is welcome to make it compatible with the upper constraint on the dark matter 
relic density. Similarly, the elastic scattering cross section for neutralinos with a 
large Higgsino component is large and can be in conflict with current experimental 
bounds. In our case, the recent constraints from XENON100 are very severe and ex- 
clude a large portion of the models we analysed. The theoretical predictions for the 
spin- independent neutralino-nucleon scattering cross-section in both BP1 and BP2 
are below the current exclusion line which, for a mass in the range 60 - 100 GeV, 
is of the order of a SI < 10~ 9 pb. Notice in this sense that increasing the singlino 
component through the decrease, e.g., of the k parameter may be useful to avoid 
this problem but this alternative might be questionable since it leads to a reduction 
in the di-photon production from the Higgs decay. 

4 Collider analysis 

Let us finally present the collider studies for the 4£ + fa signal of the benchmarks 
shown in Tables [1] and |2j As discussed in Sec.|2j the benchmark points can be 
classified in terms of the Higgs masses, namely, 

• Scenario I: (m h o,m fc o, m a o) = (126, 98, 3.6) GeV and 
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• Scenario II: (m h o,m h o, m a o) = (126, 62, 7.6) GeV. 

In these benchmark points, BR(/i!] — > a?a°) is comparable to or larger than BR(/i2 — > 
ZZ*), and BR(/i ( j ) — > a°a5) ~ 1. The light pseudoscalar boson a? decays dominantly 
into the di-tau final state with branching fraction of ~ 79 (92) % in Scenario I (II). 
Given that the pseudoscalar mass is larger than 2m M , the di-muon decay mode is 
negligible in Scenario II, while it remains sub-leading for Scenario I. Thus, restating 
Eq. (J4J), the leading final state will be 

hi, 2 -> a?a? -> r+r-r+r- -> 2£+ + 2T + # T , (5) 

with I = e, /i, or r^. Here, electrons and muons are coming from leptonic tau decays, 
while r/j denotes the tau-jet originated from hadronic tau decays. Occasional muons 
can come from a\ — > fi + n~ decay process, which is sub-leading or negligible in the 
chosen benchmarks. The source of the missing energy is associated with neutrinos 
resulting from tau decays. In order to determine the feasibility of the signal at the 
LHC run with the 8 TeV center-of-mass energy, y/s, which has recently finished 
its operation with success, we consider inclusive search channels. In other words, 
our analysis differs from past studies which usually consider one or two exclusive 
channels such as 4/x or 2/i + 2-^ with the 14 TeV beam condition. 

We have generated Monte Carlo (MC) event samples of the Higgs signals for 
a proton-proton collision at yt~s = 8 TeV using HERWIG++ 2.6.1 [83] with the 
CTEQ6L1 parton distribution functions (PDF) [SI]. h\ ~^ ZZ* — > 4£ processes are 
also generated since the final states are similar as in Eq. ([5]). The generated event 
samples have been scaled to 25 fb" 1 of integrated luminosity, which corresponds 
to the LHC data accumulated in the year 2012. To simplify the analysis, we here 
consider leading-order cross sections for all the processes. The production cross sec- 
tion and the decay width of the SM Higgs boson are calculated with higlu [85] . 
then the NMSSM Higgs cross sections have been obtained by corrections accord- 
ing to the total decay widths and BR(/i5 2 ~~ * 99) calculated with nmssmtools. 
The generator-level events are further processed with the fast detector simulation 
program Delphes 2.0.3 [SB] using a modified CMS detector card. In the detector 
simulation, jets are reconstructed by using the anti-/c 4 algorithm [87] with the radius 
parameter of 0.4, and the fe-tagging efficiency is set to be 60%. Isolated muons are 
required to have the transverse momentum > 5 GeV within the pseudo-rapidity 
range |?7| < 2.4. 

The dominant SM backgrounds include the Drell-Yan (DY), bb, cc, semi- or di- 
leptonically decaying ti, WW/Z, and ZZ/j processes. Direct Jjijj or T productions 
can in principle contribute, but they are found to be almost negligible. To estimate 
the backgrounds, bb and cc event samples are generated by Pythia 6.4 [88] and 
CTEQ6L1 PDF sets in the various pr bins, and the remaining dominant background 
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processes are generated by HerwigH — h at the matrix-element level. The parton 
showering and the hadronization are performed by HerwigH — h The MC samples 
have been fed into Delphes using the same detector card as that used for the Higgs 
signals. Then, in order for suppressing backgrounds while keeping the Higgs signals 
of interest as much as possible, the basic event selection cuts are imposed as given 
below. 

• At least two pairs of oppositely-charged leptons including the tau-jet with 
p T > 6 GeV, 

• no 6-tagged jet, 

• missing energy fa > 23 GeV. 

When selecting the leptons, the priority is given to the hard electrons and isolated 
muons since the reconstruction efficiency of the tau-jet is relatively poor@ For the 
events that passed the basic selection cuts, the combinatorial ambiguity of pairing 
the leptons is resolved by calculating all the possible M T2 defined as in Ref. [89] . 

Mt2 = niin 

k «+4 2 )= pC T 

where (i = 1, 2) are the transverse masses calculated by the measured transverse 
momenta of the charged leptons and the invisible momenta kjp of the neutrino 
system, which are decided by a numerical minimization over all possible splittings 
that maintain the missing energy condition constructed in the event. Basically the 
Mt2 has been introduced for a decay topology like a pair production of the squark 
that decays into a quark and the LSP. The decay topology given in expression (jSJ) 
is essentially different from that because there are at least four neutrinos in each 
event. However, since the pseudoscalar bosons are emitted back-to-back in the rest 
frame of the CP-even Higgs boson and the their decay products are nearly in the 
colinear direction, the neutrinos sharing the same parent pseudoscalar boson can be 
considered as one invisible particle, so that the M T2 is applicable. Since two pairs of 
oppositely-charged leptons are initially selected, there are two possible ways of the 
pairing. We choose one pairing that gives the smaller value of Mt 2 - This method 
can be justified by the fact that the M T2 value would be bounded from above by the 
parent particle mass m a o if one chose the right pairing, whereas there is no such a 
restriction if the pairing was wrong [90] ■ Since we here consider a light pseudoscalar 
boson, an upper cut on the M T2 < 10 GeV is further employed. The M T2 distribution 
is shown in Fig. [T] for both benchmark points. We can appreciate a very clear peak 
and the endpoint structure around the pseudoscalar mass for both scenarios. After 
choosing the pairing of the leptons, the transverse momentum p^ e = |Pt~ + Pt I 

4 The branching ratio of the leptonic tau decay mode is about 34%. 
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Figure 1: The Mp2 distributions of four visible leptons + missing energy signals for (left 
panel) the BP1 and (right panel) the BP2. All the event selection cuts described in the 
text have been imposed. 

and the angular separation AR,£+£- = a/ (Arje+e-) 2 + (A<p£+£-) 2 are calculated for 
each lepton pair. Since it is likely that the leptons are fairly boosted in the nearly- 
colinear direction, the selection cuts on the p 1 ^ 1 and ARg+£- can be applied as 
below. 

• p^ e > 10 GeV for each pair and 

• min{Ai?«_, Ai^_} < 0.8. 

On top of that, the invariant masses of the oppositely-charged lepton pairs are 
require to be m e + e ~ < 10 GeV to ensure that the lepton pairs come from the light 
pseudoscalar. Since there is no source of hard jets except the initial state radiation, 
the event including a jet with p^ > 25 GeV is discarded. Finally, we calculate the 
cluster transverse mass defined as 

M\ = Um? v + |p¥l 2 + ^t) 2 - |p¥ + W 2 , (7) 

where V denotes the four-lepton system [91]. The endpoint position of the Mp dis- 
tribution corresponds to the parent particle mass of the pseudoscalar bosons, i.e., 
m h o or m h o if there exists an event with vanishing invariant mass of the neutrinos 
and all the particle tracks are on the transverse plane. However, since one cannot 
distinguish the decay event of h\ from that of h\ as they share the same final states, 
and due to the limited statistics in the 8 TeV LHC experiment Jf] it is hardly see 

5 One might expect a secondary edge at m h o below the primary edge at m h a for high luminosity 
data and a fairly good reconstruction efficiency. 
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selection cuts 


BP1 


BP2 


DY 


7 7 

00 


ft 


Basic cuts 


1.29 


1.11 


3363.47 


5239.82 


272.46 


M T2 < 10 GeV 


0.64 


0.34 


2.25 


1015.97 


4.80 


p^ e , min AR1+1-, and m £ + £ - 


0.64 


0.34 


1.88 


1005.20 


4.64 


Jet-veto 


0.31 


0.19 


0.38 


8.27 


0.03 


Mt < 170 GeV 


0.17 


0.17 


0.38 


0.0 


0.03 



Table 5: Cut flow of the signals and backgrounds in fb. See the text for the detailed 
descriptions of the cuts applied. 

the edge structure near m h o in the Mt distribution. Furthermore, our MC study 
found that the efficiency of the transverse mass variable can be highly dependent 
on the reconstruction efficiency of the objects, in particular, the missing transverse 
momentum from the neutrinos. Another possible collider variable is the invariant 
mass reconstructed by the Mp2-assisted on-shell approximation of the invisible mo- 
menta [92] , which has been found to be useful for reconstructing the resonance mass 
peak in the similar decay topology [93J . However, because of the limited set of signal 
data, it is inevitable that one postpones to use the variable until the future 14 TeV 
LHC experiment with somewhat high luminosity. Still, the Mt can be used to de- 
fine the signal region. Conservatively, we set Mt < 170 GeV ~ m t . This definition 
can be quite stringent if there is a new physics process with invisible particles in 
the final state mediated by a heavy new particle which can be a background to the 
Higgs signal. 

After imposing all the event selection cuts, cc and the di-boson processes turn out 
to be completely negligible. For the estimation of the signal significance, we show 
how the cross sections of the Higgs signals and the backgrounds change under each 
event selection cut in Table [5j The cut efficiency of each scenario is rather different. 
One reason seems to be the dissimilar efficiency of the object reconstructions. In 
BP1 the decay products of the pseudoscalar boson are relatively adjacent as it is 
likely that they are more boosted than in the case of BP2, and this can affect the 
correct identification of the final-state particles. 

Once the event has been passed all the selection cuts, we construct the invariant 
mass of all visible leptons, m 2 £+2£-- The peak position of the invariant mass dis- 
tribution is generically different from m h o and m h o due to invisible neutrinos from 
the tau decays. However, one can still estimate the energy scale of the resonances 
as can be seen in Fig.|2J Although it seems to be a bit difficult to see a golden 5cr 
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Figure 2: Invariant mass distributions of four visible leptons of the signals for (left panel) 
the BP1 and (right panel) the BP2. All the event selection cuts described in the text have 
been imposed. 

discovery of the Higgs signal at the 8 TeV LHC experiment, we expect that the Mr-2 
and the invariant mass with the dedicated selection cuts described here might help 
to glimpse the non-standard Higgs scenarios proposed in this paper. 



5 Conclusions 



In the light of the recent detection at the LHC of a scalar boson with a mass of 
126 GeV and the possible hint for an excess in LEP that could correspond to a 
second, lighter scalar we investigate the phenomenology and detectability of viable 
models in the NMSSM. We consider two scenarios that contain two light Higgses, 
in which the mass of the lighter scalar is either in the range of 96 - 100 GeV or 
lighter than half of the SM-like Higgs mass. Both of them are characterized by the 
presence of a very light pseudoscalar boson, 2m T < m ai < 2mi,, a range which is yet 
to be explored further at the LHC. The Higgs phenomenology of these scenarios is 
extremely rich, involving cascade Higgs decays that lead to four leptons in the final 
state, with or without missing energy. We perform a reconstruction of these channels 
for two representative benchmark points, showing that the pseudoscalar mass can 
be successfully determined from the M T2 distribution of four visible leptons and 
missing energy. We also investigate the reconstruction of the scalar Higgs masses 
from the invariant mass of the four visible leptons. Although the presence of invisible 
neutrinos from tau decays leads to a broadening of the distributions and a shift in 
the peaks of the corresponding mass distributions, one can still estimate the scale 
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of the Higgs masses. We encourage experimentalists to search for inclusive leptonic 
decay modes in order to shed some light on these possible two Higgs scenarios. 
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